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ABSTRACT 

We have measured Johnson BV and Cousins RI magnitudes for comparison stars near five TeV 
blazars. We compare our values with published values, spanning 25 years in some cases, to identify 
those stars that are most likely proven stable. To avoid zero-point offsets mimicking long-term vari- 
ability, we based our analysis on the standard deviation between measurements after a mean offset 
between data sets was removed. We found most stars to be stable at the 0.04 magnitude level. We 
confirm two stars as variable and identify two others as possibly being variable. In each of the five 
fields there are at least two stars, and typically many more, that show no evidence of variability. 
Subject headings: AGN - comparison stars - optical photometry 



1. INTRODUCTION 

The standard model of Active Galactic Nuclei (AGN) 
successfully explains blazars as high-energy jets directed 
toward us at a small angle to the line of sight. This orien- 
tation Doppler-enhances the emission intensity, contracts 
the variability time scales, and leads to apparent super- 
luminal motion in ejected radio knots. In extreme cases 
we observe intensified synchrotron radiation from the ra- 
dio to the UVX-ray band and inverse-Compton radiation 
at higher energies, up to the TeV dom ain. 

As pointed out bv lAbdo et al.l ()2011t ) the details of the 
physical processes underlying blazar emission are still 
unknown. In their study of Mrk 421 both a hadronic 
model based on the Sync hrotron-Proton Blazar model of 
iMiicke et al.l ()200lL |2003| ) and a leptonic model based on 
the 1-zone Synchro tron Self-Compton model described in 
iFinke et al.l (pOOl were able to reproduce the observed 
spectral energy distribution, the less-constrained leptonic 
model agreeing marginally better with the data. A dis- 
similar distribution of matter in these models leads to dif- 
ferent predictions regarding the speed with which emis- 
sion intensity can vary. So a key to discriminating be- 
tween these two scenarios is the rate of multi- wavelength 
variability. Optical variability studies therefore provide 
data that help constrain these competing models. 

It is essential in optical variability studies to first estab- 
lish well-calibrated secondary standards within the field 
of view (FOV). These standards should be tied back to 
known systems such as Johnson-Cousins UBVRI to en- 
able interobject comparisons. They must be known to 
be stable so that their own intrinsic variability is not as- 
cribed to the blazar itself, especially in cases where there 
are only one or two appropriately bright standards in the 
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field. Establishing these standards now enables current 
and future all-sky surveys such as the Large-Scale Syn- 
optic Telescope (LSST; Ivesic et al. 2008) to confidently 
tie their measurements to historic data. 

Well-measured standards exist for many of the brighter 
blazars, but the set is not complete. Finder charts with 
photometric standards for 92 prominent blazars compiled 
by the Whole Earth Blazar Telescope (WEBT) project 
are maintained at the University of Heidelberg. In these 
data the epochs of most recent calibration span 20 years 
with observations often in only two or three filters. In the 
set as a whole, photometric indices exist for 19% in [/, 
57% in B, 60% in V, 78% in R and 55% in /. Most missing 
indicies are for objects that are rarely, if ever, monitored. 
Of greater concern here is that few of the standards, in- 
cluding those of the commonly observed objects, have 
been verified since the mid to late 1990s. It is important 
for monitoring accuracy to periodically verify standard 
magnitudes, weed out slowly varying stars, and extend 
the calibrations to unmeasured filters as needed. 

In their excell ent study of AGN and blazar compar- 
ison s tandards, iGonzalez-Perez. Kidger. fc Martm-LuisI 
(|2001f) make the case for establishing many secondary 
standards per field since larger numbers reduce the sta- 
tistical error. This is correct but only with regard to 
random errors; it does not improve the systematic error 
if all stars in the same field possess the same zero-point 
offset. Additionally, the overall error can increase if too 
many of the secondary stars are at large distances from 
the object of interest where flat-fielding errors tend to 
be greatest. On the other hand, if at least one star is 
proven to be unvarying and tied accurately back to a 
standard system, then the precision of the photometry 
can be enhanced from comparing the object of interest 
to an ensemble of stars while the accuracy can be deter- 

^ http: / / www.lsw.uni-heidelberg.de / projects / extragalactic/charts / 
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mined by referencing the zero-point to the weU-known 
star. 

In this paper we examine the variabihty of stan- 
dard stars surrounding five of the more significant TeV 
blazars: Mrk 421, H1426-f428, Mrk 501, 1ES1959-I-650, 
and BL Lac. We describe our observations in Section 2. 
In Section 3 we present our results and in Section 4 we 
discuss the conclusions. 

2. OBSERVATIONS 

The objects o bserved were taken fro m the northern 
blazars listed by iHoran fc WeekesI (|2004l ) as having been 
detected at TeV energies. Each field has between four 
and sixteen previously published calibrated BVRI stan- 
dards in a 15' square FOV. 

Observations were made with the BYU Robotic Ob- 
servat ory for V ariable Object Research (ROVOR) tele- 
scope (iMoodv e t al. 2012). This telescope consists of an 
RC Optical 0.40 m Ritchey-Chretien optical tube on a 
Paramount ME German-equatorial mount. The CCD 
camera is an FLI Proline PL003 with a 1024 x 1024 24 
/im pixel SITe detector. The FOV is 23.4' on a side with 
a resolution of 1.37"/pixel. The filter set was Johnson- 
Cousins BVRI c losely matched to the specifications de- 
fined bv Bessell (pTffl . 

Data were taken on thirteen photometric nights be- 
tween June fc September 2009. Primary standards from 
ILandoltl (12001) were observed up to 20 times each night 
and never less than five times. Approximately one-third 
of the time at the telescope was spent observing the stan- 
dards and a unique transformation was determined for 
each night. 

All of our nightly obse rving routines followed closely 
the procedures outhned in lLandoltl ()2007[ ). We began and 
ended each night by observing a Landolt field containing 
ten or so standard stars. We observed each standard 
field first near the meridian at an airmass no higher than 
1.2 and every hour thereafter until reaching an airmass of 
2.0. At that time we began the process over with another 
standard field near the meridian. The blazar standards 
were always observed within the airmass range covered 
by the standards. We observed only in the western half 
of the sky with the telescope on the eastern side of the 
pier to avoid flat-fielding inconsistencies from the image 
inverting as the telescope crossed the pier. The CCD 
temperature was always at -30° C to eliminate errors in 
dark current calibrations caused by non-linearities with 
temperature. We observed in a "palindrome" sequence 
of VBRIIRBV to place our mean observations in each 
filter at the same average air mass for a given setting. 
All of our images were 300 seconds in 5, 180 seconds 
V, 90-120 seconds in R and 30-60 seconds in /. Finder 
charts from images taken with the BYU West Mountain 
Observatory (WMO) 36" telescope are given in Fig. [T] 

2.1. Data Analysis 

We used two software programs to analyze our data: 
iraf and Mira. For about three-fourths of the nights 
iraf was used to apply the dark, bias, and flat-field 
frames as well as perform the aperture photometry. Mira 
was used for the remainder of the nights. To ensure that 
both programs produced similar results, two nights were 
reduced and photometered using both programs. The re- 
sulting instrumental magnitudes agreed to within 0.001 



mag. In either case we used an aperture radius of 5" and 
the stellar FWHM was typically on the order of 3". 

Our observations were place d on a standar d footing 
following the steps outlined by ILandoltl (j2007t ) . We dis- 
covered that the atmospheric extinction varied as a func- 
tion of azimuth as well as altitude on many of our nights, 
an effect especially apparent in the / band. This is be- 
cause our observatory is near a desert and on a valley 
floor, making it more susceptible to errors from nightly 
extinction variation. To be cautious, we applied a unique 
atmospheric extinction correction to each program and 
standard field on each night. Since our program fields 
were always observed across an airmass of one to two, we 
could measure the extinction for each field through that 
range. We did this by first forming a magnitude versus 
airmass relation for each comparison star in a given field 
for each night. We then found BVRI extinction coeffi- 
cients for each star in the field using a least-squares fit 
to this relation. The average of these coefficients was 
then used to remove the effects of extinction from all of 
the stars in that field. This process was repeated for all 
fields. 

The extinction-corrected magnitudes were transformed 
onto the standard system using the equations 

V = V0+ ^i,{B ~V)+Cv, 
B -V = ^ibv{b - v)o + (bv, 

V - R = ^yr{v - r)o -I- Cvr, 

V ~ I ^ liyi{v - i)o + C„. 

The /X coefficients are the standard color terms, the ( co- 
efficients are the zero points, the ca pital letters are the 
literature values from ILandoltl (|2009f ). and the lower-case 
"sub zero" terms are the instrumental magnitudes, cor- 
rected for extinction. Values for the coefficients were 
calculated using a least-squares fit to the extinction- 
corrected instrumental magnitudes and literature values 
of the Landolt standards. 

3. RESULTS 

Magnitudes of the comparison stars, associated errors 
(tr) , and number of nights the star was observed are pre- 
sented in Table [T] In this table the names of blazars 
appear in the first column, while star numbers are in the 
second column. The next four columns give magnitudes 
and errors in the listed filters, while the final column lists 
the number of nights the star was observed. The errors 
quoted are the standard deviation of our final calibrated 
values. 

We estimated stellar variability by first finding the 
magnitude difference between our measurements and 
previously published measurements for each star. To 
minimize the effects of zero-point offsets, we found the 
average difference in each filter between our data and the 
data in a given publication and subtracted it off. We then 
averaged the absolute value of the residual of each star 
across all publications, which we plot against apparent 
V magnitude in Fig. [2] These plotted points therefore 
are a normalized variation between publications. 

The error of our measurements is on the order of 0.02 or 
0.03 magnitudes which, again, includes the uncertainty 
of referencing back to the Landolt system. This is sim- 
ilar to the errors quoted in most of the previous studies 
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against which we are comparing. Adding typical errors in 
quadrature gives a value of about 0.04 magnitudes which 
is a reasonable approximation of the expected value of an 
average point in Fig. [2] Stars with residuals above this 
value, shown as a horizontal line in that figure, are sus- 
pect for variability. One can see five stars that are sus- 
pect by this criterion. In the sections below we discuss 
these outliers individually. 

The statistic formed above is better for finding variabil- 
ity than for determining stability since it has a mean off- 
set removed as explained. To establish the stars of great- 
est stability, we felt it best to simply choose those with 
the closest agreement in published photometry across all 
references with no offset removed. However, we did not 
consider references with an offset larger than 0.04 mag- 
nitudes. 

Stars of greatest estimated stability in each field are 
listed in Table[2] The names of blazars appear in the first 
column, while star numbers are in the second column. 
In columns three through six we give their BVRI mag- 
nitudes and errors as determined by an error-weighted 
mean across all references. The errors thus calculated 
are from random sampling only, no attempt was made 
to formally estimate the syst ematic errors . We fol- 
low the numberi ng schemes of iVillata et al.l ()1998[ ) and 
iFiorucci fc Tostil 11996) in designating the comparison 
stars. Other identifications are listed in column 7 for 
clarity. 

3.1. MRK 421 

The star field around MRK 421 is sparse and 
dominated by the 6th magnitude star 51 UMa 
and its 12th magnitude companion BD+39 2414B. 
iGonzalez-Perez. Kidger. fc Martm-LuisI (|2001l ) observed 
Mrk421 but there were no other stars inside their ap- 
proximately 7' square FOV. There are, however, five 
published comp arison stars outside this distan ce from 
IFiorucci fc ^ti (1996) Doroshenko et al.l (|2007[ ). 

Star 3 is problematic. It is the biggest outlier in Fig. 
O Our uncorrected BVR values for it are 0.28, 0.04, 
and 0.15 ma gnitud es brighter than those presented by 
IVillata et al.l (pM) . This may be an intrinsic variability 
or just an effect from the orientation of our telescope 
secondary mirror spider arms, which, as can be seen in 
the finder chart, cause difi'raction spikes from 51 UMa to 
cut through it. Either way, we recommend that star 3 
not be used as a comparison star. 

Excluding star 3, our values averaged across all fil- 
ter bands agree to within 0.035 magnitudes with the 
average values of the other four stars presented by 
IFiorucci fc Tostil (119961 ) and within 0.006 magnitudes 
with Doroshenko et alT (|2007l ) . Since this agreement is 
within the published errors, we conclude that stars 1,2, 
4, and 5 show no evidence of variability. 

3.2. H1426+428 

The field near 111426+428 is also sparsely populated. 
Because of this, we have adde d only one n e w com parison 
star to the four presented bv lSmith et al.l (|1991| ). 

With only one set of published values to compare 
against, this field cannot be examined in the same man- 
ner as the others, so we simply compare data. When 
averaging differences across all bands, the iSmith et al.l 



(fl99l values differ from ours by -0.017 for #1, 0.0625 
for #2, 0.038 for #3, and 0.13 for #4. The differences for 
stars #1, and #3, are within the errors or nearly so. Star 
7^3 shows a more significant difference but we cannot as- 
cribe it to variability from this one data point. However, 
star ^4 is so different in value tha t it is hard to imag ine 
how errors of observation by either [Smith et all (|1991fl or 
us could cause such a discrepancy. Even when removing 
the average difference between our two papers, the offset 
is still 0.07 magnitudes, the value we chose to plot in Fig. 
121 Upon that discrepancy alone we recommend avoiding 
star 4. 

3.3. MRK 501 

Mrk 501 has a nicely populated field with eight pub- 
lished magnitud es fro m three sources Fiorucci & Tos^ 
(1996), Villat a et al.l (119980. a nd Doroshenko et al] 
(2005). .Gonzalez-Perez. Kidger. fc Martfn-Lu is {Mm 
publish magnitudes for 18 faint stars encompassing the 
faintest ones near the object. We publish values for three 
more. 

Stars 1 and 4 show excellent agreement across all au- 
thors. The difference between our values and the sources 
cited above are 0.033, 0.016, -0.004, and -0.004 for star 
1 and -0.01, 0.01, 0.008, and -0.009 for star 4 respec- 
tively with t h e cav eat that the B value for star 4 from 
IVillata et all (|1998[ ) was rejected as an obvious outlier. 
Being close to Mrk501, they are excellent references, es- 
pecially the brighter star 1. 

For star 2 the values from IFiorucci fc Tostil (|1996D are 
internall y consistent but o f f by 0.07 magnitudes from 
ours. In iDoroshenko et al.l (|2005D they are again inter- 
nally consistent but disagree with our values by 0.02 mag- 
nitu des. Star 3 shows the same trend; our values disagree 
with IFiorucci fc Tosti ( 1996.) bv 0.11 magnitude s and by 
0.04 magnitudes with IDoroshenko et al.l (|2005l ). These 
differences could be a zero-point offset or they might 
be intrinsic var iation, particularly with star 3 since our 
agreement with IDoroshenko et all ()2005D data is usually 
much better than 0.04 magnitudes. Since the accuracy 
of these stars are in question, and stars 1 and 4 are so 
solid, we recommend referencing to stars 1 and 4. 

3.4. 1ES1959+650 

1ES 1959-K650 has a well-populated field. IVillata et all 
(11998") pu blish B and V values for seven stars and 
Doroshen ko et al.l (j2007[ ) publish BVRI values for the 
same set. We add four more stars. 

Star 5 is actually the variable MM Draconis. Its vari- 
ability is apparent both in Fig. [5] and in the large inter- 
nal error of our data (se e Ta b le lit. We also confirm the 
finding of Doroshenko et al.l (j2007D that star 3 is vari- 
able and should not be used. All remaining stars show 
no obvious variation. The mean difference between our 
values and those of stars 1, 2, 4, 6, ai id 7 have remark- 
ably g ood agreement with residuals of IDoroshenko et all 
(|2007 !). These residuals averaged over all filters are 0.011, 
0.007, 0.005, .001, and .013 respectively. The good agree- 
ment is probably from the large number of observations 
in their study (33) and ours (10). 

3.5. BL Lac 

The field of BL Lac is also well populated and well 
calibrated, with a variety of stars close to the blazar. 
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We consi der here the c ompar is on star magnitudes pub - 
lish ed by iSmith et alj (ITOSl. iFiomcci fc Tostil (Il99l) , 
and lGonzalez-Perez. Kidger. fc Marti'n-Luisl (|2001|V 
All standardized stars in this field show 

no obvious variation. We c oncur with 

iGonzalez-Perez. Kidger fc Martin-LuisI (I2001D that 



the photoelectric values from ISmith et all (I1985D are 
systematically brighter and exclude them from our 
analysis here and in Table [2] The errors averaged across 
all filters are on the order of 0.02 magnitudes for all 
stars with star 5 having the largest variation of 0.035. 

4. CONCLUSION 

Most stars examined show no evidence of variability 
to within the errors of observation. Further, there are 
within each field at least two demonstrably stable stan- 
dards stars, suitable for referencing. That these are often 
the brighter, closer stars unsurprisingly illustrates that 
minimizing error from photon statistics and flat-fielding 
produces the best data. 



Referencing against many stars provides the most 
precise photometry. But for many applications a ref- 
erence against the known stable stars will be suffi- 
cient and simpler. We note that Doroshenko et al. 
(2005, 2007) tied their photometry back to the Johnson- 
Cousins system through what they considered to be the 
most reliable measurements from other papers, pa rticu- 
larlv lGonzalez-Perez. Kidger. fc Martfn-LuisI (|200H ). We 
found that overall our numbers agreed best with theirs, 
often to within than 0.01 magnitudes. 
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15'x 15' •• . BL.Lac •. 

Fig. 1. — Finder charts for the blazars with comparison stars identified. North is to the top and east is to the left. All charts are 15' x 
15' as indicated. All images were taken with the WMO 36" telescope in the V band. 
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Fig. 2. — A plot of magnitude difference residuals against V magnitude for all comparison stars considered in this paper. The larger 
the residual, the more likely the star is variable. The stars well above the nominal random error line of 0.04 magnitudes are Mrk 421 #3 
(0.122), H1426+428 #4 (0.079), and 1ES1959+650 #3 (0.075) and #5 (0.056). 
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TABLE 1 

BVRI MAGNITUDES OF COMPARISON STARS. 



Blazar 


Star 


B(<t) 


V(^) 


R(a) 


I(^) 


n 



MRK 421 



14.97 
16.15 
16.41 
15.17 
14.44 



0.02) 
0.03) 
0.05) 
0.02) 
0.03) 



14.34 
15.56 
15.73 
14.13 
13.58 



0.03) 
0.04) 
0.03) 
0.07) 
0.05) 



0.02) 
0.05) 
0.06) 
0.01) 
0.01) 



0.02) 
0.02) 
0.02) 
0.05) 
0.03) 



13.97 (0.03) 
15.19 (0.05) 
15.09 (0.06) 
13.52 (0.03) 
13.06 (0.03) 



13.63 
14.82 
14.61 
13.05 
12.65 



13.26 (0.03) 
14.12 (0.03) 
12.97 (0.02) 
15.05 (0.04) 
14.74 (0.02) 



H 1426+428 



15.60 
15.39 
14.19 
16.00 
15.98 



14.15 
14.55 
13.39 
15.39 
15.18 



12.49 
13.73 
12.60 
14.72 
14.36 



0.02) 
0.02) 
0.05) 
0.04) 
0.06) 
0.10) 
0.07) 
0.02) 
0.03) 
0.02) 
0.02) 



0.02) 
0.02) 
0.03) 
0.03) 
0.05) 
0.04) 
0.04) 
0.02) 
0.02) 
0.02) 
0.02) 



12.08 (0.02) 
12.70 (0.02) 
14.69 (0.03) 
14.92 (0.02) 

14.97 (0.03) 
14.91 (0.03) 
14.61 (0.02) 
12.55 (0.02) 

13.98 (0.02) 
13.03 (0.02) 
11.30 (0.02) 



MRK 501 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 



13.54 
14.03 
15.87 
15.93 
16.14 
16.77 
16.71 
13.48 
14.94 
14.37 
12.81 



12.60 
13.17 
15.13 
15.31 
15.41 
15.64 
15.44 
12.90 
14.34 
13.54 
11.82 



11.63 
12.28 
14.30 
14.57 
14.56 
14.36 
13.94 
12.19 
13.63 
12.56 
10.83 



0.02) 
0.03) 
0.03) 
0.03) 
0.14) 
0.04) 
0.04) 
0.03) 
0.02) 
0.05) 
0.02) 



0.02) 
0.01) 
0.02) 
0.01) 
0.13) 
0.02) 
0.02) 
0.01) 
0.01) 
0.02) 
0.02) 



12.28 (0.01) 
12.55 (0.01) 
12.33 (0.02) 
14.04 (0.02) 
14.00 (0.10) 
14.77 (0.03) 
14.75 (0.03) 
11.88 (0.01) 
13.10 (0.01) 
14.72(0.03) 
13.93 (0.02) 



1ES1959+650 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 



13.34 
13.43 
15.02 
15.26 
15.65 
15.94 
15.96 
12.69 
14.07 
16.19 
15.14 



12.69 
12.90 
13.30 
14.50 
14.64 
15.21 
15.22 
12.21 
13.49 
15.27 
14.38 



11.90 
12.22 
11.47 
13.61 
13.43 
14.38 
14.30 
11.57 
12.73 
14.27 
13.53 



0.03) 
0.03) 
0.04) 
0.06) 
0.03) 
0.03) 
0.03) 
0.03) 
0.03) 



0.02) 
0.03) 
0.03) 
0.03) 
0.03) 
0.03) 
0.03) 
0.03) 
0.03) 



11.98 (0.03) 
13.77(0.02) 
13.72 (0.03) 
15.06 (0.03) 
13.84 (0.02) 
12.71 (0.03) 
12.84 (0.03) 
13.77 (0.03) 
13.61 (0.03) 



BL Lac 



14.65 
15.14 
15.74 
16.33 
15.48 
14.31 
13.99 
15.05 
15.07 



12.94 
14.26 
14.45 
15.54 
14.46 
13.31 
13.28 
14.25 
14.16 



11.15 
13.31 
13.07 
14.58 
13.30 
12.17 
12.41 
13.26 
13.07 
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TABLE 2 

Estimation of photometric values for. demonstrably stable stars. 



Field 


Star 


B{a) 


V(a) 


R(<7) 


I(<7) 


Other Designations 




1 


14.982(0.009) 


14.384(0.006) 


14.022(0.004) 


13.692(0.003) 








2 


16.173(0.010) 


15.571(0.007) 


15.200(0.004) 


14.866(0.005) 


35 




MRK 421 


4 


15.135(0.010) 


14.124(0.006) 


13.534(0.004) 


13.018(0.004) 








5 


14.390(0.012) 


13.571(0.006) 


13.062(0.004) 


12.628(0.004) 


25 






1 


15.61(0.02) 


14.16(0.01) 


13.25(0.01) 


12.46(0.01) 






H 1426+428 


2 


15.45(0.02) 


14.60(0.01) 


14.15(0.01) 


13.76(0.01) 


32 






3 


14.20(0.01) 


13.40(0.01) 


12.99(0.01) 


12.61(0.02) 








1 


13.540(0.003) 


12.598(0.003) 


12.083(0.003) 


11.613(0.001) 







MRK 501 


4 


15.961(0.018) 


15.309(0.010) 


14.931(0.009) 


14.572(0.010) 









1 


13.359(0.012) 


12.686(0.012) 


12.280(0.004) 


11.916(0.004) 








2 


13.446(0.014) 


12.888(0.008) 


12.534(0.004) 


12.217(0.004) 






1ES1959+650 


4 


15.277(0.014) 


14.501(0.008) 


14.038(0.004) 


13.619(0.004) 








6 


15.968(0.015) 


15.204(0.012) 


14.758(0.004) 


14.365(0.004) 








7 


16.001(0.015) 


15.225(0.012) 


14.738(0.004) 


14.315(0.004) 








1 


14.643(0.020) 


12.936(0.011) 


11.966(0.014) 


11.105(0.016) 




12^^ 




2 


15.178(0.012) 


14.278(0.008) 


13.760(0.011) 


13.313(0.013) 




2lie 


BL Lac 


3 


15.723(0.017) 


14.452(0.004) 


13.695(0.013) 


13.024(0.009) 




1616 




4 


16.380(0.021) 


15.549(0.009) 


15.051(0.009) 


14.561(0.010) 




I7I6 




5 


15.484(0.026) 


14.439(0.007) 


13.797(0.010) 


13.237(0.013) 


3I6 






6 


14.327(0.023) 


13.298(0.004) 


12.689(0.006) 


12.161(0.015) 


25I6 






8 


15.097(0.008) 


14.254(0.016) 


13.755(0.017) 


13.262(0.020) 


28'' 






9 


15.144(0.013) 


14.191(0.020) 


13.590(0.017) 


13.070(0.023) 


29^ 





Note. — (1 ) [Poroshenko ot aL (2005), (2) Smith ct aL (1991), (,3) Villata et aTl ( 11993 '). ('4') ISmith eTal] ( 119851) . (5) 
IBertaud et al.l ( Il96ai , fei , Gonzalez-Perez. Kidgcr. &: Martm-Luia (,200l1) 



